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We demonstrate that a suspended nanowire forming a weak link between two superconductors can 
be cooled to its motional ground state by a supercurrent flow. The predicted cooling mechanism has 
its origins in magnetic field induced inelastic tunneling of the macroscopic superconducting phase 
associated with the junction. Furthermore, we show the voltage-drop over the junction is propor- 
tional to the average population of the vibrational modes in the stationary regime, a phenomena 
which can be used to probe the level of cooling. 
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Nanoelectromechanical systems (NEMS) are fast ap- 
proaching the limits set by quantum mechanics [iMSl- 
Achieving such conditions requires that the mechanical 
subsystem can be brought into, and detected, in its quan- 
tum mechanical ground state. In general this condition 
demands that an energy quanta associated with the me- 
chanical motion is much larger than the energy associated 
with the thermal environment. For an oscillator with a 
mechanical frequency of 100 MHz this implies tempera- 
tures as low as a few mK. However, using oscillators with 
higher mechanical frequencies the quantum limit can be 
reached, as recently demonstrated by O'Connell et al. [4J . 

The most common device geometries of NEMS to date 
consist of mechanical oscillators in the form of cantilevers, 
suspended beams or microtoroids. These typically have 
much lower resonance frequencies than those reported in 
Ref. [4], hence reaching the quantum limit in these de- 
vices is very challenging. To circumvent this problem, 
back-action cooling of the mechanically compliant ele- 
ment is often employed whereby the number of mechan- 
ical vibrons is reduced without necessarily lowering the 
ambient temperature. Suggestions for different cooling 
mechanisms are plentiful, see e.g. Refs. [E|-[9|. Common 
to these is that the oscillator is cooled either by coupling 
its mechanical oscillations to electromagnetic photons or 
a flow of charge carriers. 

In the present paper we suggest a new mechanism of 
cooling not previously considered and show that ground- 
state cooling of the mechanical oscillator is possible. Con- 
sidering the nanomechanical oscillator as a weak link in 
a current-biased Josephson junction we show that we 
can access a regime analogous to the resolved side-band 
limit [10], whereby the number of mechanical vibrons in 
the system can be reduced by a factor of ~100. In the 
limit of a high mechanical quality factor the resulting vi- 
bron population is shown to be well within the quantum 
regime. 

The cooling mechanism considered here is achieved by 
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FIG. 1. (Color online) Schematic diagram of the system. 
(Left) A suspended nanowire of length L forms a weak 
link between two current-biased superconducting leads. The 
transverse magnetic field H is applied perpendicular to the 
nanowire. (Right) The equivalent electronic circuit. A con- 
stant current I is applied to the Josephson junction which is 
connected in parallel to a capacitor C and a resistance R. 



coupling the mechanical vibrations of the oscillator to 
the supercurrent through the junction. Below we show 
that the suggested setup not only allows for ground-state 
cooling of the mechanical oscillator, but simultaneously 
probes the macroscopic nature of the superconducting 
phase associated with the junction. As such, the pro- 
posed system allows for interesting physical observations 
on both the mechanical and the electronic subsystems. 

Figure Q] shows a schematic picture of the system con- 
sidered. It consists of a metallic carbon nanotube sus- 
pended over two superconducting leads biased at a cur- 
rent /. Transverse to the in-plane motion of the nanotube 
a magnetic field H is applied which induces coupling be- 
tween the bending modes of the wire to the supercurrent 
through it [ll|. Below we analyze the influence of the 
electromechanical coupling and show that for resonant 
current-biased conditions this may lead to ground-state 
cooling of the vibrations of the nanowire. 

In our analysis we restrict the description of the me- 
chanical degrees of freedom of the nanowire to the funda- 
mental bending mode, which is considered as a harmonic 
oscillator with frequency uj. The Hamiltonian describing 
the system presented in Fig. [1] has the form, 

H = AE c n 2 - jhcj) - Ej cos(</> - <M) + hwb ] b . (1) 
Here, n is the operator for the number of Cooper pairs 



on the junction and <\> is the corresponding operator for 
the superconducting phase [l2| h] = i). In ([1]), 
E c — e 2 /(2C) is the Coulomb energy where C is the 
capacitance of the junction, j = I/(2e) is the flow of 
Cooper pairs and Ej is the Josephson energy. The op- 
erators [b] are creation [annihilation] operators for 
the oscillator where u — b + W is the dimensionless 
deflection of the wire. In the above, the parameter 
$ = 4g7TLHu zp /^o characterizes the strength of cou- 
pling between the mechanical and electronic degrees of 
freedom. Here, u zp = (H/(2muj)) 1 ^ 2 is the zero-point 
amplitude of the nanowire, m and L is the effective mass 
and length of the suspended part of the wire respectively, 
$o = 7rh/e is the flux quantum and g is a numerical fac- 
tor of the order of unity which accounts for the profile of 
the fundamental bending mode [l3j . 

The third term in (pQ) describes on the one hand the 
Lorentz force on the nanowire induced by the Joseph- 
son current. On the other hand, it gives the deflection- 
dependence of the Josephson current due to the motion 
of the wire in the magnetic field In what follows 

we consider a nanotube of length L ~ 1/im, for which 
u zp < 1A, in a magnetic field H ~ IT. With these pa- 
rameters $ < 0.3, and we consider only the linear terms 
in the expansion of (pQ) with respect to With this 
expansion the Hamiltonian reads, 

H — Hj + H 7 n + Hint > (2) 

Hj = 4:E c n 2 — jhcj) — Ej cos cp , H m = fruPb , 
Hint = -Ej${fr + b) sin0. 

Here, Hj is the Josephson Hamiltonian, which under 
the condition j < Ej/h, describes the electronic sub- 
system in the so-called tilted washboard potential. In 
(j2j), Hint describes the interaction between the mechani- 
cal and electronic subsystems with H m the Hamiltonian 
of the former. 

Below we will take the Columb energy to be much 
smaller than the Josephson energy, AE C /Ej <C 1. This 
condition implies that the characteristic interlevel dis- 
tance between the quantized states of the Josephson 
junction associated with a given local minimum of the 
washboard potential, Huj p = (SEjEc) 1 ^ 2 , is much smaller 
than the height of the barrier separating different local 
minima. Here, uj p is the plasma frequency of the junc- 
tion. We also take the external temperature T to be 
low, T < hup/hs, such that transitions between states 
associated with different local minima can only occur 
through under-barrier tunneling. A schematic diagram of 
the quantum state of the electronic subsystem described 
through Hj is shown in Fig. [2j 

Under-barrier tunneling between two consecutive val- 
leys in Fig. [2] changes the state of the Josephson junction 
through the associated change of the phase. Such tunnel- 
ing events, commonly referred to as macroscopic quan- 
tum tunneling (MQT), are greatly enhanced if the two 
energy levels involved in the transition are in resonance. 
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FIG. 2. (Color online) Schematic diagram of the tilted wash- 
board potential U(<f>) = —Ej cos cp — jhcj) as a function of phase 
cp at current-bias / = e/it(uj p — uj). Here, / labels the valleys of 
the potential and a =f, ^ are the two energy levels within the 
valleys considered. In the above, <E>T is the inelastic tunneling 
amplitude between two energy levels in consecutive valleys. 
The quantity T is the transition rate from the second to the 
first level within a valley generated by interactions with the 
quasiparticle environment (see text). 



This can be achieved by tuning the current-bias. Thus, 
we define the critical bias current /* as the current which 
ensures that the lowest (first) level in a given valley is res- 
onant with the second level in the next valley, /* ~ ecu p / it 
[l4j . As the potential defined by Hj is only to first ap- 
proximation parabolic, the spacing between the energy 
levels within a given valley is not constant. As such, we 
will in the following only consider tunneling between the 
two lowest electronic states and neglect any coupling to 
higher levels. This is justified as the, e.g., the second 
and third levels are far from resonance if the junction is 
biased at I ~ /* (see Fig. EJ 

The electronic system in Fig. [2] is coupled to the me- 
chanical subsystem by the magnetic field. As such, 
MQT can in the present situation also be accompa- 
nied with the emission/absorption of a quanta of me- 
chanical energy, huj. Performing a WKB analysis for 
the MQT amplitude we find that the overlap integrals 
for the inelastic channels is of the order of where 
T cx Hujp exp(— tt(Ej /(2E C ) 1 ^ 2 ) < Huj is the tunneling 
amplitude in the elastic channel. Here, we note that the 
^-dependence of Hint only leads to a renormalization of 
the parameter g in the definition of Also note that due 
to the large separation in energy, uj <C uo pi the electrome- 
chanical coupling will not introduce additional tunneling 
channels between the higher electronic energy levels. 

The inelastic tunneling channels change the number of 
mechanical vibrons such that cooling of the oscillator is 
possible if transitions through the absorption channel can 
be promoted. Below we show that this can be achieved 
by tuning the bias current so that the absorption channel 
is resonant; the first level in a valley I is separated by fvjj 
from the second level in Z + l as shown in Fig. [21 A further 
condition for cooling is that the electronic subsystem, 
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once in the second energy level, relaxes to the lower level 
at a rate T which is faster than the rate at which the 
system tunnels back with the emission of a vibron, T > 
T/h. Such relaxation arises due to interaction with the 
quasiparticle environment as discussed further below. 

To perform a quantitative analysis of the system we 
introduce the basis |Z, a) where I labels the valleys of the 
potential and a =t, I labels the energy levels inside a 
given valley (| is the first and f is the second level). In 
this basis the Hamiltonian reads, 

n = n + n r , (3) 

H =nj + nm=Yl + tuJtb)% a) (I, a\ , 

l ,n,cr 

U r = T + gt ) + l) I* + 1> t) (I, I I + hx. . 

i 

In the above, = fux) p ma- — Inhl/e are the eigenvalues 
for the electronic degrees of freedom in the basis \l,a), 
where = 1 and m± = 0. From the form of the Hamil- 
tonian (|3]) one can see that due to the electromechanical 
coupling the number of vibrons in the system is not con- 
served and may change due to macroscopic tunneling of 
the electronic system from one valley to the next. 

To describe the joint dynamics of the electronic and 
mechanical degrees of freedom we will start our analysis 
from the Liouville-von Neumann equation for the density 
matrix p of the system, 



dp 
&t 



%o +Ht,P 



J(f>Y 



7(1 + n B )C h (p) + jn B C hi (p) . 



(4) 



Here, J(p) is a phenomenological damping operator for 
the electronic system [14j , 

J (fi = ~l (E l*'t>(i,t \p + p\IA){IA + 
r^2\i,i)(iA\p\i'A)(i',i\- 



(5) 
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In the equivalent circuit scheme (see Fig. [T]) this damp- 
ing derives from the parallel resistance R, which in the 
present situation causes the system to decay from the t 
state to the 4- state in a given valley. In (J5j), T = uo p /Q e i 
is the electronic damping rate, where Q e \ — u p RC is the 
corresponding quality factor. Here we consider Q e i ^> 1 
which implies that the influence from the electronic quasi- 
particle environment on the tunneling processes is neg- 
ligible [14l4l6j. We will further suppose that the quality 
factor Q e i is so large that broadening of the second energy 
level, Aujp = uo p / (2Q e i), is small enough for the inelastic 
resonance transitions to be resolved, Auj p < uj. 

The second damping term in (J2J), £>a(p) = (2apa) — 
a)ap — pa^a)/2, is the standard Lindblad operator which 
models interactions between the oscillator and the ther- 
mal environment. Here, 7 = uo/Q is the mechani- 
cal damping rate with Q the quality factor and n B = 



(ex.p(f3hu) — where /3 = (/c^T) -1 , is the average 

number of vibrons in thermal equilibrium. 

Below we investigate the stationary solution to (j4]). To 
find this solution we perform a standard perturbative 
analysis in the small parameters T/(hT) : j/T oc e <C 1 
and look for a solution of the density matrix of the 
form p = po + epi + e 2 p2... (for a full derivation of 
the results presented below see Appendix [XJ). Substi- 
tuting this into (j4]) one finds that the leading order solu- 
tion po has the form p = Y^i, n K> h n )Po(M, n)(hU n \, 
where the index n labels the Fock state of the oscilla- 
tor. From (J2J) we also find the first order correction 

Pi = E/,n,j=-l,0,l \ l + 1 ^ Tl + 3) C 3& n )( l ^M + h ' C - 

where the sum j^/Cj(Z,n) = Cj(n) satisfy the following 
relation, 



Cj(n) 



T} n) P(n) T(n) 
t/ST/2' j 
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J 
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Here, P(n) = Po(h >U n ) 18 the population of the vibra- 
tional modes of the oscillator. Developing the perturba- 
tive expansion one finds that the equation for the second 
order term, ^2, can only be resolved if P(n) satisfy the 
following equation, 

(r_+7(l + n B ))[(n + l)P(n + 1) - nP(n)} + 

(r+ + jn B )[nP(n - 1) - (n + l)P(n)} = . (6) 

Here, Tj are the different tunneling rates; j = — , 0, + are 
respectively the absorption, elastic and emission channel, 



4^27-2 



4T 2 



4(Aj-±) 2 + /i 2 r 2 ' 

AT = ^7+i,t — • 



A T± = J-q ± Hlj . 



Considering the operator for the potential over the 
Josephson junction V = i\H,4>]/(2e) (in our represen- 
tation (j> = 2tt G |Z, a) I (I, cr|) we find, 

^ = J E 7 " + ^) + !. t>(i,4- I + hx. . (7) 

This implies that the stationary bias voltage, V = 
Tr(Vp), is zero to leading order in p. Thus, the potential 
drop is given by the first order correction to the density 
matrix, V = Tr(Vpi), which implicitly depends on the 
coefficients Cj(n). Solving equation ((6]) we find that the 
average number of vibrons, (n) = J2 n n P( n ), is given by, 



(n) 



7 + r_ - r+ 

and that the voltage drop scales with (n) as, 
nh 



V 



(r_(n) + r + r + ((n) + i)) . 



(8) 



(9) 



Here we note that the potential drop in the stationary 
regime is primarily determined by the elastic tunneling 
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FIG. 3. (Color online) Average vibron population (solid) and 
bias voltage (dashed) in the stationary regime as a function 
of the current bias. Here, $ = 0.3, r = cj/4, T = hu/20, 
tie 20 and Q =10 5 . 



Appendix A: Derivation of density matrix 

The evolution of the density matrix is governed by the 
Liouville-von Neumann equation (J2J), 



dp 

at 



7 (1 



+ J(p)+ 
>yn B £ ht (p) . 



(Al) 



In what follows we will consider the stationary solution 
of (|Aip by performing a perturbative analysis in the 
small parameters T/(hT)^j/T <C 1. In particular we 
will consider the limit of high mechanical quality factor 
Q such that 7 = u/Q < T/(ti). To start the analy- 
sis we take the total density matrix to be of the form, 
P — Po + £pi + e 2 p2... and equate powers of e. With this 
we find the following equations, 
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0(e°) (A2) 



rate, To. This is consistent with the physical processes 
discussed, i.e. in the limit 7, T + — >> we get (n) = 
(complete ground state cooling as no heating channel is 
open) and V oc r (the system moves down the tilted 
washboard potential at the rate r which conserves the 
number of vibrons). 

In Fig. [3] we plot both the average stationary popu- 
lation of the mechanical subsystem and the correspond- 
ing voltage drop as a function of the bias current. As 
expected, the lowest occupation is achieved when I = 
I* — eu/7r (see Fig. [2]). In this regime, we find that ground 
state cooling of the mechanical subsystem is possible if 
the resolved side-band limit, uj > T, is achieved. Under 
conditions when the bias current is I > I* the tunnel- 
ing events discussed above will lead to pumping of the 
mechanical subsystem, in which case the above analysis 
does not apply once the limit T((n) + 1) ~ HT is reached. 
This regime will be discussed in future work. 

To conclude we have shown that a suspended nanowire 
which forms a weak link in current-biased Josephson 
junction can be cooled to its motional ground state. This 
effect derives from the coupling of the mechanical motion 
of the nanowire to the electronic degrees of freedom by 
a magnetic field. Furthermore, we have shown that by 
operating the system under optimal bias-current condi- 
tions the occupation factor of the vibrational modes can 
be greatly decreased. Also, we have found that the po- 
tential drop over the junction might be a sensitive probe 
of the stationary vibron population as it scales with the 
average number of vibrons. 

This work was supported in part by the Swedish VR 
and SSF and by the EC project QNEMS (FP7-ICT- 
233952). 
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, 0(e x ) (A3) 



0(e 2 ) (A4) 



Solving the above equations at each order of e we find 
po = n \l,l,n)po(l,l,n)(l,l,n\ which satisfies (jA2|) . 
Similarly, the first order correction to the stationary den- 
sity matrix is determined from (jA3|) as, 

pi = ^2 I' + l?tj n + i) c j(l, n )(l,l, n \ + n - c - 1 



-1,0,1 



Cj(l,n) = 



7- (n) Po(U^) 
-AJ 7 j + ihT/2 ' 



(A5) 



Substituting this into ()A4|) we find the equation for the 
coefficients po by tracing out the spin (t,i) degrees of 
freedom, 



- oi 4(A^-) 2 + ^r 2 



Tj n - j) \ 2 po(l-l,Un-j)- 



|7I (n) | 2 p (U,n; 



j(l + n B )[npo(l,l, n) - (n + l)p (l,U n + 1)]+ 
jn B [(n+ i)po(l,l,n) - np (/,|,n- 1)] . (A6) 

Tracing out the valley index I we recover the expressions 
presented in the paper, i.e. equation (jA5|) gives 



-AFj +ihT/2 
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whereas equation (jA6|) gives, 

(r_+7(l + n B ))[(n + l)P{n + 1) - nP{n)] + 

(r+ + 7 %)K(^ - 1) - (n + l)P(n)] = . (A7) 



In this expression the relationship between the coeffi- 
cients are, 



r± = r 



4(AJ* ± ) 2 + h 2 T 2 ' 



r n = r 



4T 2 



4(A7b) 2 + h 2 T 2 ' 



^ (n) (QH-C*(n))=P(n)r j iV, 




N 



In the above we note that (|A7|) gives the balanced equa- 
tion for the probability P(n) of finding the oscillating 
nanowire in the state n. The stationary average distribu- 
tion of the vibrational modes is then given by the solution 
to this equation, 



(n) = nP(n) 



n=0 



r_ + 7 -i\, 



The density matrix pi allow us to evaluate the potential 
drop over the junction in the stationary regime. Follow- 
ing the derivation outlined in the paper we find that the 
lowest order term of the density matrix, po, does not con- 
tribute to the potential drop as it is diagonal in the spin 
basis. As such, the potential drop is uniquely determined 
from pi. 
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